Interplay of Kondo and superconducting correlations in the nonequilibrium Andreev 

transport through a quantum dot 
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Using the modified perturbation theory, we theoretically study the nonequilibrium Andreev trans- 
port through a quantum dot coupled to normal and superconducting leads (N-QD-S), which is 
strongly influenced by the Kondo and superconducting correlations. From the numerical calcula- 
tion, we find that the renormalized couplings between the leads and the dot in the equilibrium 
states characterize the peak formation in the nonequilibrium differential conductance. In particular, 
in the Kondo regime, the enhancement of the Andreev transport via a Kondo resonance occurs in 
the differential conductance at a finite bias voltage, leading to an anomalous peak whose position 
is given by the renormalized parameters. In addition to the peak, we show that the energy levels 
of the Andreev bound states give rise to other peaks in the differential conductance in the strongly 
correlated N-QD-S system. All these features of the nonequilibrium transport are consistent with 
those in the recent experimental results [R. S. Deacon et ai, Phys. Rev. Lett. 104, 076805 (2010); 
Phys. Rev. B 81, 12308 (2010)]. We also find that the interplay of the Kondo and superconducting 
correlations induces an intriguing pinning effect of the Andreev resonances to the Fermi level and 
its counter position. 

PACS numbers: 73.63.Kv, 74.45. -l-c, 72.15.Qm, 73.23.-b 



I. INTRODUCTION 

Electron transport through nanofabrications has at- 
tracted much attention in the studies of fundamental 
quantum physics as well as potential future devices. In 
particular, a quantum dot (QD), which has discrete en- 
ergy levels where the electrons are correlated, provides an 
ideal arena to study the local Coulomb interaction effect 
on the transport^"— . The magnetic doublet states with 
spin 1/2 are stabilized at an isolated QD with an odd 
number of electrons and the strong Coulomb interaction, 
which results in the Coulomb blockade for the transport 
through the QD coupled to leads. At sufficiently low 
temperatures, however, the local moment of the doublet 
states is screened by the electrons of the leads owing to 
the Kondo effect, and thus the Kondo singlet is stabi- 
lized, resulting in an anomalous enhancement of the zero 
bias conductance^"—. 

If we replace the leads by s-wave superconductors, a 
different situation arises; the doublet is not screened due 
to the lack of low-lying energy states of the leads. Even 
in the system, a singlet state can be stabilized due to 
the superconducting proximity effect, and the system 
thus shows a transition between the doublet and the sin- 
gletiSr— . Away from the transition point, one of the two 
states becomes the ground state and the other an ex- 
cited state which is localized at the QD, i.e., the Andreev 
bound state. In this system, however, it is difficult to di- 
rectly observe the Andreev bound states via transport 
measurements because of a supercurrent and a multiple 
Andreev reflection process^"— . 

Recently, Deacon et. al. have observed the Andreev 
bound states experimentally not in the above system 
but in the system with a QD coupled to normal and 



superconducting leads (N-QD-S) where an Andreev re- 
flection dominates the transport^i^. In the N-QD-S 
system, however, the doublet states should be replaced 
by the Kondo singlet state owing to the screening by 
the electrons of the normal lead (N-lead), leading to a 
crossover between the Kondo singlet to the supercon- 
ducting singlet. A lot of studies have thus far focused 
on how the competition between the Kondo and super- 
conducting correlations affects the Andreev transport ex- 
perimentally22r— and theoretically^"—. Indeed, Kondo- 
type anomalous phenomena have been observed in the 
measurement of zero bias conductance in the recent ex- 
periment^. 

Experimentally, characteristics of the Andreev bound 
states emerge under nonequilibrium steady-state condi- 
tions where a finite bias voltage is applied to the N-lead. 
Some theoretical studies have dealt with the nonequilib- 
rium transport properties in an N-QD-S system with em- 
phasis on the influence of the Kondo effec1j^"-2^'^"i^'^"ii 
and also on the Andreev bound states^. However, the 
coexistence of the phenomena related to the Kondo ef- 
fect and the Andreev bound states in the experiments 
indicates the necessity of further theoretical studies; it is 
needed for the comprehensive understanding of the trans- 
port to include the Andreev bound state as well as the 
interplay between the Kondo and superconducting corre- 
lations into the theory. 

In this paper, we study the nonequilibrium Andreev 
transport, by taking into account the above different as- 
pects of the N-QD-S system in a unified way. To this end, 
we employ the modified second order perturbation theory 
(MPT) used previously by Cuevas et al.— and extend it 
to the nonequilibrium steady-state conditions. The MPT 
was originally formulated in the equilibrium Anderson 



model^, then has been used in several different systems, 
e.g., a quantum dot coupled to normal leads^"— and as a 
impurity solver for the dynamical mean-field theory^^^. 
Furthermore, we exploit the exact solution of the QD- 
S system with an infinitely large superconducting gap, 
which still has the essence of the Andreev bound states, 
to improve the perturbation theory. By systematically 
examining the nonequilibrium transport properties in a 
wide variety of the system parameters, we demonstrate 
that the theoretical results obtained in this paper are 
qualitatively in agreement of the recent experimentsi^i^ 
We note that a part of the present results was briefly re- 
ported in ref. 1471 . 

This paper is organized as follows. In Sec. II, the 
model Hamiltonian is introduced and we formulate the 
modified second-order perturbation theory in Keldysh- 
Nambu space of the Green's function. Section III, we 
assess the validity of our method in the equilibrium case 
and define the renormalized parameters which character- 
ize the electron transport in the nonequilibrium states. 
The results of nonequilibrium transport are shown in Sec. 
IV. We also analyze the superconducting pair amplitude 
and the local density of states at the QD in the nonequi- 
librium states. The correspondence between the theo- 
retical and experimental results is also discussed in this 
section. A summary is given in Sec. V. 



II. MODEL AND METHOD 

A. Model Hamiltonian and Keldysh Green's 
function in Nambu space 

In order to describe the electron transport in the N- 
QD-S system, we use a single level QD coupled to a nor- 
mal metal and a superconductor, which is applicable for 
the system with large level spacing of the QD, 

H = i?QD + ^^N + ^^S + -f^TN + ^^TSi (1) 

where 

-ffqo Ed ^ "do- + C/ Ud-fUdi, (2) 

^fN = ^(e^-/^iv)cLc,.., (3) 

ka 

qa q 

i?TN = ^(iNcLrf.+H.C.), (5) 

ka 

J^TS = E(^s4o4 + H.c.). (6) 
go- 
Here, creates an electron with spin a at the QD which 
has an energy level €d and the Coulomb interaction U. 
Here, Uda = d^da- cl.^ (ajo-) denotes the creation oper- 
ator of an electron with spin a and wave vector k (q) in 



the normal (superconducting) lead. The superconduct- 
ing lead is assumed to be described by the BCS Hamilto- 
nian with a superconducting gap As — Aexp{i9s)- The 
QD is coupled to the normal and superconducting leads 
labeled by a = A^, S* with hybridization ta ■ 

In order to define nonequilibrium steady states of 
the N-QD-S system with a bias voltage V, we treat 
the Coulomb interaction J7 as a perturbation. The 
non-interacting problem can be solved exactly with the 
Keldysh Green's function technique in Nambu space, 
from which we can define the chemical potentials of the 
normal and superconducting leads as /ijv = and 
/is = 0, respectively. 

In the noninteracting case, several different Green's 
functions in Nambu space at the QD are defined as 

nUf ^[^tW'4(^')])o {[d4t).d^{t')]}o 

^ ^ ^ '\{[di{t),di{n])o {[diit),d,{t')])o 
^ \{[dl{tld\{t')]), {[dl{t),d,{t')]),) 

n<(f f')^ ( i4it')d^{t))o {d^{nd^{t))o\ 

^ > \{diit')di{t))o (d,(t')di(i))or 



'{d^{t)dl{t'))o {d^{t)d^{t'))o^ 
{dl{t)dl{t'))o {dl{t)d^{t'))o^ 



( 



where and g"^ denote the retarded and advanced 
Green's functions, which are also used in the equilibrium 
case, and g< and g^ represent the lesser and greater 
Green's functions. We consider a sufficiently wide band 
of electrons in the leads, in which the coupling strength 

rAr(s)(w) = MtN{S)\T,k{q) ^(^ - '^k{q?^ becomes a con- 
stant ^N(S)- Integrating out the electron degrees of free- 
dom in the two leads, we obtain the Fourier transformed 
Green's functions, 

g-{u) = {{u + i7^)I - EdtTs - Sr(c^))"' , (11) 

g^iij) = {{lj - trj)I - ed(T3 - S?(a.))-^ , (12) 

g<{uj) = -g^{i,)^<ioj)g^{i,), (13) 

g>{cj)^-g^{Lo)i:>{oj)g^{u), (14) 



where 



-*(rAr + rs/3(^)) *rs/3(w)^ 



(16) 



f{uj - flN) 







f{uj + /ijv) 



-z2rsRe[/3(c^)] ( J^'VH' (17) 
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butions to the retarded and advanced self-energies, 



-'2ndV 



FIG. 1: Second order self-energy diagrams we consider in this 
paper. The solid line indicates the propagator and the dashed 
line the Coulomb interaction. 



Here, i] is a positive infinitesimal and cr^ (i — 
1,2,3) is a Pauli matrix in Nambu space. = 

If we obtain the self-energies due to the Coulomb inter- 
action, the full retarded and advanced Green's functions 
are determined from the Dyson equation. 



(19) 



The full lesser and greater ones are calculated from the 
relation, 

G<->{u;) ^ -G''(w)[s<'>(a;) -|-S<'>(w)jG°(a;). (20) 

Here we examine the self-energy Sfj, using the pertur- 
bation theory in the Keldysh Green's function formal- 
ism. The first order contributions to the retarded and 
advanced self-energies are S^j," = U{N) with 



(21) 



where (n^) denotes the expectation value of the electron 
number at the QD per spin. There is no first-order con- 
tribution to the lesser and greater self-energies because 
the Coulomb interaction at the QD takes place without 
delay. The second order skeleton diagrams are depicted 
in Fig. [TJ The corresponding contributions to the lesser 
and greater self-energies are obtained from the equations, 

S2<dH - ~U'I ^n<(.. + c.i)T2 {g>{u^^))'^cT,i22) 
where 

n<>(a.) = j ^(^g<^>ico,)g<,>{co-u,) 

-9f^>{ui)9^{^{uj-oj,)), (24) 

where gij denotes the (i,j) component of g. Using the 
above self-energies, we calculate the second order contri- 



27r 



duJi 



[S2<dK) 



sr„d(^i)] 



(25) 

Although the second order self-energies are believed 
to give reasonable results for the nonequilibrium trans- 
port through a strongly interacting QD coupled to two 
normal leads (N-QD-N system) at least for the particle- 
hole symmetric cas o^^'^^ , this technique is not directly 
applicable to the N-QD-S system because of the lack of 
symmetry of the leads. Main difficulty in our N-QD-S 
system comes from the fact that the simple second or- 
der self-energies do not correctly give the formula in the 
"atomic limit" where the QD and the leads are discon- 
nected. Indeed, the qualitatively correct description of 
the Kondo effect in a particle-hole symmetric N-QD-N 
system with the second order self-energy is ensured by 
the fact that the corresponding formula becomes exact 
not only in the weak-C/ but also in the atomic limit In 
the N-QD-S system, furthermore, the superconducting 
correlations at the QD, which come from the supercon- 
ducting proximity effects, must be taken into account, so 
that we have to introduce a suitable "atomic limit" to 
study the strong-C/ regime. For this purpose, we here 
make use of the exact solution of the QD-S system with 
an infinitely large superconducting gap. 

Below, we describe how to construct the modified self- 
energies, S2nd, which reproduce the correct results in 
the atomic limit as well as the weak-U limit within the 
above second-order perturbation. The self-energy due to 
the Coulomb interaction up to second order thus reads 

S[/ = Slst + S2nd- 



B. Superconducting atomic limit 

In the limit of A — )■ oo, the quasiparticle degree of 
freedom in the superconducting lead is decoupled from 
the QD. Therefore, the Hamiltonian is simplified as. 



^ - ^dot 



Hn + Htn 



H.c. 



(26) 

Und^;ndi{2^) 



where = T s Gyi]i{i6 s) ■ In this case, Ts corresponds 
to the effective superconducting gap at the QD owing to 
the proximity effect. 

In the limit of Tn — > 0, the QD is decoupled from the 
normal lead and the effective Hamiltonian (pS)) becomes a 
one-site problem with the superconducting paring poten- 
tial and the Coulomb interaction. Hereafter, we call the 
limit of (A — >■ oo,Tm/U — 0) "superconducting atomic 
limit"—. In the superconducting atomic limit, the self- 
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energies at the QD can be exactly obtained as, 

i::U^)^U\{{u + zv)I-A)-\ (28) 
i:t,^iu:)^U\{{uj-zrj)I-Ar\ (29) 

S<tm(^) = -^'^VU'X {{UJ + ^V)I - A)-' 
x{{I-<T,)/2 + {N)) 

X {{u ~ ir^)I - Ar\ (30) 
S>^(c.) = z277;72x ((c. + irj)I - A)-' 
x{{I + cTs)/2-{N)) 
X {{uj - irj)I - A)-\ (31) 

where 

X ^ {n,){l-{n,))-\{d^d^)\\ (32) 



A*^-U{d^d^)* -Ed - (7(1 - (nd)) 



(33) 



Note that we omit the first order contributions of U in 
the atomic-limit self-energies. 

Next, we consider the second order self-energies, fol- 
lowing the formula derived in the previous section. We 
assume that one-particle Green's functions in the second 
order diagrams (Fig. [T]) are dressed with energy shifts, 
which are determined by the following one-body Hamil- 
tonian, 

H = iJdot + Hn + Hs + Htn + Hts, (34) 
i?dot = (ed + C/(n^))^d^d^ 



+ {U{did-^)dldl+H.c.), 



(35) 



where {Wd) and (did-^) are the effective parameters repre- 
senting the energy shifts of the one-particle Green's func- 
tion. 

In the superconducting atomic limit, the above one- 
particle Green's functions behave like a S function, so 
that the second order self-energies can be evaluated as, 

SLdH^C^'xo((c^ + *r;)/- Ao)-\ (36) 

^^2nA^)^U\o{{u^-^v)I-^o^\ (37) 

^2ndi^) ^ -^^vU\o ((^ + iv)I - Ao)"' 
x((I-rT3)/2+{N)o) 

X {{lj - ir^)I - Aor\ (38) 

5^2''nd(^) ^ i27jU\o ((^ + iv)I - Aq)"' 
x((J + ^3)/2-(Ar)o) 

X ((a;-^77)7- Ao)"\ (39) 



where 



Xo = {nd)o{l-{nd)o)-\{did^)o\^ (40) 

- yA*, + u{d^d^y -ed-u{m))' ^^'^ 



/^r\ = f Mo {did^)o\ 



(42) 



Here, {nd)o and (d^(i^)o are the expectation values of 
the particle number per spin and the superconducting 
correlation at the QD under the one-body Hamiltonian 

(El. 

We find that in the superconducting atomic limit, the 
second-order self-energies have the functional forms simi- 
lar to the exact ones, except that they have different con- 
stants; (xo, Aq, {N)o) and (x, A, (AT)). Exploiting this 
fact, we construct the modified second order self-energies 
in the following. 



C. Modified second order perturbation theory 

First, we formulate the modified self-energies for the 
retarded and advanced sectors. For the sake of clarity, we 
follow the procedure of the modified perturbation theory 
(MPT) in the N-QD-N system by Kajueter and Kothar-^. 
In this procedure, the modified self-energies are assumed 
to have the following functional forms. 



SJ„,(a.)=A[[SJ„,(a.)]-^-B 



(43) 



where A and B should be determined for the self-energies 
to reproduce the exact ones in the high-energy limit as 
well as the superconducting atomic limit. The exact 
high-energy limit of the self-energies can be calculated 
from the continued-fraction expansion of the correspond- 
ing Green's function a^^'^" . 



s[;"(c^) = [/(Ar) + ^/- 



0(A) 



(44) 



The first term coincides with the first order self-energy 
in U. On the other hand, in the limit of w — ^ oo, the 
modified self-energies are expanded as. 



S2„d(^) = . I 



(45) 



The coefficient A in eq. p5|) is determined from the con- 
dition that the leading terms of the modified self-energies 
are identical to the corresponding ones in the exact self- 
energies in eq. (El) . Accordingly, we set A as x/xo- 

We next determine the matrix B from the condition 
that the modified self-energies give the correct values in 
the superconducting atomic limit. In the limit, the mod- 
ified retarded and advanced self-energies become 



KM ^ U'xii^ ± - Ao - U\oB)- 



(46) 



In order to eliminate the difference between the r.h.s of 
the above equation and that of eqs. ([25)1 and (|29l) . we set 
B as follows, 



B 



1 



1 - {nd) - {nd) -{did^) - {did^)\ 
UxQ K~{did^)* - {did^y -1 + (nd) + {nt} J 



(47) 
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Note that in the hmit of A — or Fg — 0, the su- 
perconducting correlations at the QD vanishes and the 
ofF-diagonal terms in eq. (|43l) become zero. Further- 
more, A = (ud) (1 - (nd)) /((n-d)o (1 - ('^d)o)) and the 
non-diagonal terms of B vanish. As a result, the modi- 
fied self-energies of eq.(|43|) just coincide with those in the 
previous studies^S. Therefore, we believe that the mod- 
ified self-energies obtained here are proper extensions of 
those used in the N-QD-N system. 

We have to calculate the modified lesser and greater 
self-energies in order to obtain the transport properties. 
By generalizing the strategy used previouslj^, we define 
the modified lesser self-energy as. 



X in." 



>nd 



(48) 



We now confirm that the above self-energy reproduces 
the atomic-limit form. Multiplying the modified self- 
energy (|48)) on the left and right by the inverse matrices 

of S2nd('*^) ^2nd('*^)' ^^^^ supercouductiug 
atomic limit. 



^^((/-<T3)/2+(Ar)o). 



(49) 



In the limit, the above matrix does not depend on uj. In 
a similar way, we multiply the atomic-limit lesser self- 
energy (|48l) on the left and right by the same matrices. 
The resulting matrix also becomes a constant in the limit. 



SLd(^) 



S<tm(w) SU(C.) 



((/-^3)/2 + (iV)). 



(50) 



Therefore, the difference between eq. and eq. ^ 
can be ignored in the limit of — ^ 0, and eq. pS)) repro- 
duces the atomic-limit form indeed. 



S2<d(^) 



,{io) (A^oo,Fjv^O). (51) 



We also define the modified greater self-energy as, 
S>,(^) = ls^„d(^)[S2'-„dH]"'S>,(a.) 



(52) 



which gives an appropriate form in the limit. 

So far, we have formulated the modified retarded, ad- 
vanced, lesser and greater self-energies. However, these 
four self-energies are not independent, but have to satisfy 
the following equality, 

S< ,(c.) - S> ,(c.) = S^„,(c.) - S-„,(c.). (53) 

We indeed confirm that the modified self-energies satisfy 
the equality. 



D. Current conservation and consistency of the 
energy shifts 



Using the resulting full Green's function, we calculate 
the current though the N-QD-S system. The current 
fiowing in the normal (N) and superconducting (S) leads 
can be calculated from the time evolution of the particle 

number operators Nn,s in each lead: /Ar(t) = 



, djVw(t) 
' dt 



and is{t) — —e — Since we assume that the system 
is in a nonequilibrium steady state, the expectation val- 
ues of these operators are time-independent, which are 
given by 



In) - 



4er 



N 



4eFc 



(54) 
-Im / dw 



Ac 



2ps{uj)fiio)Gl^{Lu) + p*{u)GU^) 



(55) 



where psii-^) = R.e[/3(cj)] and Gij denotes the com- 
ponent of G. We define the current / in this system as 
/ — {In) = — {Is)- However, it is known that the current 
calculated by the second order perturbation theory may 
not be conserved in some quantum dot systems except 
in a special condition^2i^iSi. In the N-QD-S system, the 
simple application of the second order self-energy usually 
breaks the current conservation rule, i.e. {In) + {Is) 0- 
In our modified second order perturbation theory, the 
problem in the current still exists. In order to resolve this 
difficulty within our framework, we introduce the source 
term A coupled to the current operator and add the term 
into the one-body Hamiltonian (p4)) . 



(56) 



HIn+Is)- 



The effective parameters {nd), {d^d^) and A are de- 
termined by the following consistency conditions on the 
energy shifts and the current conservation. 



U{nl) = U{nd)+Rc[-E''2ni{nN)] 

^2nd(Ais)]l 
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U{did^) = U{did^) 
{In) + {Is) = 0, 



(57) 



where [5]2jjjj]ij denotes the (i, j) component of the modi- 
fied retarded self-energy, {ud) and {did^) are also deter- 
mined in a self-consistent manner. 

Here, we check the J7 — > limit of the modified self- 
energies. In the small-[/ limit, Sjnd becomes the sim- 
ple second order self-energy Sjnd which is calculated 
from the one-particle Green's functions dressed with the 



6 



mean-field energy shift because the consistency condi- 
tions in eq ([57)1 are reduced to (n^) = {rid) = {nd)o 
and {did^) — (d^rf^) = {did^)o- Therefore, S 



-'2nd 



IS 



also reduced to Sjj^^j evaluated with using the mean-field 
Green's functions. 

Note that the above modified self-energies are applica- 
ble to impurity systems with or without superconducting 
correlation. We will show below that the above method 
works very well except for some special cases with large 
bias voltage where we cannot find the convergent parame- 
ters (fid), {d^d^} and A. In this paper, we mainly focus on 
the reasonable parameter region where the bias voltage 
is not so large. We demonstrate that a variety of intrigu- 
ing phenomena emerge due to the interplay between the 
superconducting correlation and the Kondo effect, some 
of which indeed reproduce the experimental results qual- 
itatively well. 



III. LINEAR-RESPONSE CONDUCTANCE 
AND PHASE DIAGRAM 

In this section, we study the transport properties in 
the linear-response regime and check the validity of our 
approximation for the electron transport. In addition, 
the renormalized couplings of tunneling are introduced, 
which clearly specify various regimes appearing in the 
noncquilibrium electron transport addressed in the next 
section. 



A. Zero bias conductance and the renormalized 
couplings in the equilibrium states 

Let us first consider the zero bias conductance obtained 
in two different ways within the same framework of mod- 
ified perturbation theory (MPT) to confirm the consis- 
tency of our approximation. Here, we concentrate on the 
symmetric coupling case, T^/Ts = 1, with particle-hole 
symmetry, ed/U — —0.5, in the equilibrium state. 

We first obtain the zero bias conductance by directly 
differentiating the current by the bias voltage. In this 
case, we have to calculate the lesser and greater self- 
energies in order to obtain the current from eqs. ([M)) 
and ([55]) . The current- voltage (IV) characteristics thus 
obtained for several values of U, are shown in Fig[21[a). 
In addition to the suppression of the current, we can see 
the enhanced nonlinear behavior. In order to observe the 
nonlinearity in more detail, we show the conductance, 
I/V, near the zero bias voltage in Fig[5l^b). In this fig- 
ure, the conductance curve for U/Tm = is almost flat 
near the zero bias voltage, implying that the linear re- 
sponse theory can be safely applied in this finite voltage 
region. For U/Tn = 5, the conductance is suppressed, 
yet keeps the fiat structure. With further increase in 
U , however, the conductance shows a convex curve; the 
linear response regime is restricted to the very tiny volt- 
age region, i.e., \eV\/A < 0.01 for U/Tn = 20. There- 



u/r^=o 



CD 




(b) 
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FIG. 2: (Color online) (a) Current-voltage characteristics for 
several values of U: Fs/Tn = 1, €d/U = -0.5, A/Fjv = 0.5 
and ksT/TN = 0.005. (b) Conductance as a function of the 
bias voltage. The parameters used are the same as in (a). 



fore, theoretical studies only on the zero bias conductance 
are not enough to understand the transport properties in 
the actual experiments in the strong Coulomb interaction 
regime. 

In the case oi ed/U = —0.5, we have an alternative 
expression for the zero bias conductance at absolute zero 
in terms of the renormalized couplings as^i^ 



d/ 
dF 



166-^ 



(rs/rjv)^ 



rs/r 



N 



(58) 



where Fjv and Fg are the renormalized couplings defined 
by 



Fat — zFat, 
Fs = z{Ts 



[S[;(0)]i2), 



Ts d[S[,(^)]n 
^ lAI duj 



(59) 
(60) 

(61) 



Here, denotes the (j,j) component of the Nambu 

matrix of the retarded self-energy at the QD. It is worth- 
while to note that Fat and Fg can be calculated only from 
the retarded (or advanced) self-energy in the equilibrium 
states and there is no need to calculate the lesser and 
greater self-energies. 

Figure [31(a) shows the zero bias conductance obtained 
in the above-mentioned two different ways. The values of 
the conductance obtained from the differentiation (trian- 
gle) well coincide with those obtained from Eq.([58|) with 
the renormalized parameters (circle). This fact confirms 
the consistency of our MPT treatment at least around 
the zero bias voltage. The consistency is assured by the 
effective parameter A introduced for the current conser- 
vation in the MPT framework. As pointed out in the 
previous section, a simple second order perturbation may 
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FIG. 3: (Color online) (a) Zero bias conductance as a function 
of U obtained from the two different methods for Ts/Tn ~ 1, 
ea/U = -0.5, A/Fn = 0.5, ksT/TN = 0, which are respec- 
tively denoted by circles (scheme of renormalized parameters) 
and triangles (scheme of direct differentiation), (b) Renormal- 
ized couplings Fat and Fg. The parameters used are the same 
as in (a). 



break the current conservation law; there are two differ- 
ent definitions of the current. In and Is- If the current 
conservation law is violated, at least one of the values of 
the zero bias conductance calculated from Im and Is is 
not consistent with the one obtained with the renormal- 
ized parameters. Therefore, the confirmation done here 
is important to obtain the sensible results for the current 
at finite bias voltage. 

In Fig. El^a), the zero bias conductance decreases with 
increasing U, which is due to the suppression of the An- 
dreev reflection by the Coulomb interaction. The sup- 
pression of the Andreev reflection between the QD and 
the S-lead and also the single-electron tunneling between 
the QD and the N-lead are seen in the Coulomb inter- 
action dependence of the renormalized couplings (Fig. 
[Sljb)). With increasing U, Fg decreases more rapidly 
than Fjv, indicating that the entire N-QD-S system is ap- 
proximately decoupled into two parts in the low energy 
region: S-lead and QD-N systems. Therefore, the Kondo 
singlet state becomes dominant in the ground state of 
the QD for large U, leading to the suppression of the 
Andreev reflection. 

We next discuss the zero bias conductance in the 
large A region, in comparison with the results obtained 
with the numerical renormalization group (NRG) calcu- 
lation^. Figure |4] shows the zero bias conductance as 
a function of U for several values of A. We note that 
similar calculations have been done by Cuevas et aL- - 
Let us first look at the case of infinitely large gap, where 
the closed and open circles denote the conductance ob- 
tained with MPT and NRG calculations. Although our 
approach is based on the perturbation expansion in U, 
the MPT results reproduce the NRG results in both weak 
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FIG. 4; (Color online) Zero bias conductance as a function 
of U for several values of A. The closed and open symbols 
denote the results of MPT and NRG. The other parameters 
are Fs/Fjv = 1, ed/U = -0.5 and fcsT/Fjv = 0. 



and strong U regions since in the MPT framework the 
effective parameters (n^) and {d^d^) are self-consistently 
determined to reproduce the atomic limit correctly. Only 
in the intermediate region around U/T^ — 6, we see 
some discrepancies between the two results (less than 0.3 
e^/h). 

Let us now discuss how the zero bias conductance de- 
pends on the superconducting gap A in Fig. |4l With 
decreasing A, the conductance for finite Coulomb inter- 
action is enhanced. It should be noted that the renor- 
malization factor z is determined not only by A and T s 
but also by the retarded self-energy due to the Coulomb 
interaction, as seen in eq. (|6T|) : for finite A and Fg, z is 
smaller than unity even in the noninteracting case. For 
small A, the renormalization by the Coulomb interac- 
tion is weak, leading to the enhancement of the zero bias 
conductance, as compared with the case of A = oo. 

We now look at the local density of states (LDOS) at 
the QD with and without the Coulomb interaction U . 
The LDOS in the noninteracting case is shown in Fig. 
[SJa). For Ts — 0, there is a broad resonance around 
the Fermi energy due to the decoupling of the QD from 
the S-lead, which means that the QD is in the mixed 
valence regime. For small Fg, the weight of the LDOS 
is suppressed at w = ±A since at the same energies, the 
divergence of DOS of S-lead occurs. With further increas- 
ing F5, the LDOS develops a pseudo gap at the Fermi 
energy and a double-peak structure appears inside the 
gap owing to the superconducting proximity effect; the 
superconducting (SC) singlet state becomes dominant at 
the QD. Since the two resonances inside the gap are re- 
duced to the Andreev bound states for Fjv/Fs = 0, we 
here refer to them as the Andreev resonances. The An- 
dreev resonances are located at cj ~ ±Fs with the same 
width Fat. The change from a single resonance to the 
Andreev resonances clearly characterizes the crossover in 
the dominant couplings at the QD, which occurs around 
Ts/Tn = 1. 
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FIG. 5: (Color online) (a) LDOS for several values of Fs: 
[/ = Ed = 0, A/Fjv = 0.5 and fcsT = 0. (b) LDOS for 
several values of U: Fs/Fjv = 1, ed/(7 = -0.5, A/Fjv = 0.5 
and feflT = 0. The inset is the enlarged picture in the region 
around the Fermi energy. 



In Fig. EJb), the LDOS at Fs/Fat = 1 (crossover 
regime) is shown for several choices of U . With the 
increase in J7, the Andreev resonances are merged into 
a single resonance, indicating that the superconducting 
correlations are reduced by the strong Coulomb interac- 
tion and the Kondo correlations are enhanced instead; 
the Kondo singlet state dominates the SC singlet state 
at the QD. The broad peaks corresponding to the charge 
excitations are also observed at w ~ /7/2 for large U . 
The U dependence of the LDOS in Fig. [SJb) is consis- 
tent with the preceding MPT calculations by Cuevas et 
ali^, though they did not address its relationship to the 
Kondo and SC singlet states. 

In the particlc-holc symmetric case [ed/U = —0.5) 
with symmetric couplings F5/F7V = 1, the strong 
Coulomb interaction favors the Kondo singlet ground 
state, as discussed above. If we change the ratio Fs/Fjv, 
however, the SC singlet state can be dominant in the 
ground state in the strong Coulomb interaction regime. 
Such examples are shown in Fig. [6l where the zero-bias 
conductance and the corresponding renormalized cou- 
plings are plotted as a function of Fg for Tn/U = 0.05. 
For Ts/U = rN/U = 0.05, the N-lead is strongly coupled 
to the QD; Fat > F5. As F5 increases, Fg increases more 




FIG. 6: (Color online) (a) Zero bias conductance as a function 
of Fs for Tn/U = 0.05, ea/U = -0.5, A/U = 0.025 and 
kBT/Tisr = 0. (b) The renormalized parameters as a function 
of Vs. 



rapidly than F jv , although both of them are enhanced be- 
cause the Coulomb interaction effects are suppressed by 
the superconducting proximity effects. We can indeed 
see that the crossover in the dominant couplings occurs 
around Ts/U ~ 0.1. Further increase in Fg leads to 
the enhancement of F5 and the suppression of F^v, driv- 
ing the system into the S-lead dominant coupling regime 
where the SC singlet is dominant at the QD. 

Since the zero bias conductance has the maximum 
value for Fs/Ftv = 1 (see eq. ([55|) ). the conductance 
shows a peak structure around Ts/U ~ 0.1 as shown in 
Fig. [HJa). Away from the crossover regime, the conduc- 
tance decreases both in the N-lead and S-lead dominant 
coupling regimes. For any finite values of Tn/U and 
A/[/, the crossover in the dominant couplings occurs at 
a finite Ts/U. We will see in the next section that the 
renormalized quantities Fjv and Tg also characterize the 
differential conductance even at a finite bias voltage. 



B. Phase diagram 

We summarize the results for the equilibrium N-QD- 
S system in the phase diagram specified in terms of the 
dominant couplings. Figure [7l[a) shows the phase dia- 
gram of the particle-hole symmetric N-QD-S system as 
functions of logarithms of Fjv/t/ and Ts/U . In this fig- 
ure, the solid lines, which are determined by Fjv = Tg, 
characterize the crossover behavior. In the left (right) re- 
gion divided by the crossover line, the coupling between 
the QD and the N-lead (S-lead) is dominant, F^r > F5 
{Tn < Ts). 

First, let us look at the case of A = 00 in Fig. 
[71[a). For large Tjq/U, the crossover line approaches 
r^/rs = 1 denoted by the dotted line because the ef- 
fects of the Coulomb interaction become weak there and 
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FIG. 7: (Color online) (a) Phase diagram for the particle-hole 
symmetric N-QD-S system in equilibrium conditions. The 
dominant state for finite Fjv and the ground state for Fat = 
at the QD are denoted in italic, (b) Phase diagram for the 
particle- hole symmetric QD-S system for Fiv = 0. The first 
order transition points of the ground state are denoted by the 
diamond (A/C/ = 0.025) and the square (A/C/ = oo). 



then Fjv/rs ~ Tn/Ts (see eqs. ^ 'And In the 

region of large Tn/U but small Ts/U, the QD is strongly 
coupled only to the N-lead and the Coulomb interaction 
is week, so that the QD is in a mixed-valent singlet state. 
With increasing F^, the proximity effects are enhanced 
and a pseudogap is formed in the LDOS at the QD, thus 
leading to the SC dominant state at the QD. On the 
other hand, in the small Tn/U region {ttTn/U < 1), the 
crossover line considerably deviates from the noninteract- 
ing one (Tn/Ts = 1) due to the renormalization effects 
by the Coulomb interaction. For ttTn/U < 1, the mixed 
valence state is gradually changed into the Kondo singlet 
state, so that the crossover from the Kondo singlet state 
to the SC singlet state occurs asTs/U increases for small 
Tn/U. 

The crossover line terminates at Ts/U — 0.5 in the 
limit of Tm/U — >■ as shown in Fig. [71[b), where the 
crossover is changed to a doublet-singlet transition. This 
is because the N-QD-S system is completely divided into 
two parts at Fjv — 0, the N-lead and the QD-S system. 
This doublet-singlet transition is easily seen in the case 
of A = cx), where the effective Hamiltonian of the QD-S 
system is simplified since the coupling between the Bo- 
goliubov quasiparticles and the QD vanishes; It has a 
single level with the superconducting pairing potential 
characterized by the hybridization Fg as already noted 
in eq. ([27l) . The resulting effective Hamiltonian can be 
diagonalized by the Bogoliubov transformation, leading 
to four eigenstates: two singly occupied states with spin 



1/2, I t) and | J,), and two states with total spin consist- 
ing of a linear combination of the doubly occupied and 
empty states. In the particle-hole symmetric case, the 
spin-0 singlet states are given by 

1^1) = ^(|0> - I n)), (62) 
1^2) = -L(|0) + |n)). (63) 

Note that the singly occupied states are degenerate (zero 
energy), and \S1) and 15*2) have the different energies, 
Es\ = ^^ — Ts and Es2 = ^^ + ^s■ Therefore, the candi- 
dates for the ground state are the magnetic doublet state, 
|ct), and the SC singlet state, |<S'l). Either of these two 
states can be the ground state depending on the param- 
eters, and a first order transition occurs at Ts/U — 0.5, 
which is denoted by the black square on the F jv = line 
in Fig.[7i;b); for Ts/U < 0.5 {Ts/U > 0.5), the ground 
state is the doublet state (singlet state). If Fat is small 
but has a finite value, the local moment of the doublet 
ground state is screened by the electrons in the N-lead 
and the Kondo singlet state becomes the ground state. 
Therefore, the characteristic behavior in the crossover 
line for small T^/U reflects a remnant of the doublet- 
singlet transition at F^r = 0. 

When the superconducting gap A becomes finite, the 
system is not so much simplified because of the existence 
of the coupling between the quasiparticles in the S-lead 
and the QD even for F^r = 0. Therefore, the competition 
between the Kondo effect and the superconducting prox- 
imity effect becomes important. Even in this case, there 
is still a doublet-singlet transition, which is confirmed by 
several authors in the problem of a magnetic impurity 
in superconductors^"— and the — tt transition of the 
QD-Josephson junctionaiSril. The transition point shifts 
toward lower F5 with decreasing A. We denote the tran- 
sition point for A/U = 0.025, which is obtained with the 
NRG calculation, by the diamond in Fig. El^b). 

Summarizing, the ground state of our system is always 
in the singlet phase for finite Fjv, where three different- 
type singlet regions are smoothly connected to each other 
via crossover behaviors. Only for Fat = 0, there exists a 
transition between the singlet and doublet states. 



C. Andreev bound states 

Here some comments are in order on the nature of ex- 
cited states. We start with the F^r = and A = 00 
case. Since there are only four discrete eigenstates at 
the QD in this case as discussed above, the excited 
states are localized at the QD. For Ts/U > 0.5, \a) be- 
comes the first excited state with the energy ut = \Esi\. 
In contrast, for Ts/U < 0.5, 151) becomes the ex- 
cited state with the energy ujf,. The other singlet state, 
1 5*2), is always the second excited state with the energy 
■^62 = £-52 — iniii(i?sii 0) which is larger than uji,- The one 



10 



particle excitation from the ground state to the excited 
states locahzed at the QD is observed as sharp peaks 
in the LDOS. These sharp peaks correspond to the An- 
dreev bound states. Therefore, there may be four An- 
dreev bound states in the LDOS when the ground state 
is a magnetic doublet with energy zbwf, and ±Wh2. On the 
other hand, when \Sl) becomes the ground state, there 
are only two peaks with energy ±ujb because there is no 
one particle excitation from \S1) to \S2). 

In the finite A case, there still exist the Andreev bound 
states inside the gap which correspond to |cr) or \S1) for 
any values of A though the energy of the bound states 
ujb cannot be obtained easily. Moreover, the second ex- 
cited state corresponding to \S2) may be outside of the 
gap and be absorbed into the continuum energy spectrum 
for smaU A^^. In that case, there are only two Andreev 
bound states at the QD even though the magnetic dou- 
blet state is the ground state. 

In the recent experiments^^i^, Deacon et al. have 
found that the Andreev bound states at the QD can be 
detected in the nonequilibrium transport measurements 
in an N-QD-S system. In the experiments, there are only 
two kinds of peaks which correspond to two kinds of the 
first excited bound states. This fact may be attributed 
to the small superconducting gap prepared in the exper- 
iments. This issue will be addressed in the next section 
focusing on the nonequilibrium transport. 



IV. NONEQUILIBRIUM ELECTRON 
TRANSPORT 



In this section, we study the nonequilibrium electron 
transport in the N-QD-S system with a special focus on 
the influence of the Kondo effect and the Andreev scatter- 
ing on the nonlinear transport. In particular, we clarify 
the origin of the characteristic structures in the conduc- 
tance profile in comparison with the recent experiments. 

Before elucidating the Coulomb interaction effects on 
the nonequilibrium electron transport, it is instructive to 
discuss the differential conductance, dl/dV, in the non- 
interacting case. Here, we set = 0, where dl/dV has a 
symmetric profile with respect to the ^ = axis. Figure 
[5] shows the differential conductance for several values 
of Ts /Tn as a function of the positive or negative bias 
voltage V. 

In the noninteracting case, the ratio of the bare cou- 
plings directly determines the nature of the system. We 
can see the dl/dV profiles characteristic in the N-lead 
and S-lead dominant coupling regimes in Fig. [UJa) 
and (b), respectively. In the symmetric coupling case, 
^s/^N = 1, the differential conductance has a maximum 
value Ae^/h at V = 0. The zero bias conductance is sup- 
pressed both for Ts/Tn < 1 and Ts/Tn > ^- In the 
nonlinear regime, however, the differential conductance 
in the two regimes has different characteristics. 

In the N-lead dominant coupling regime as seen in Fig. 
[SJa), the differential conductance is suppressed in whole 
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FIG. 8: (Color online) (a), (b) Bias voltage dependence of the 
differential conductance for several values of Ts, U/Tn = 0, 
ed = 0, A/Fjv = 0.5 and ksT/TN = 0. Here, we only show 
the results either for positive or negative V since dl /AV is 
symmetric with respect to = 0. (c) Peak position of dl /dV, 
Vc, in comparison with Ts- 



subgap voltages with decreasing Fg because of the sup- 
pression of the Andreev reflection. In the Fg — limit, 
the sharp peaks appear at the gap edges and the pro- 
file of dl/dV becomes similar to the density of states 
in the S-lead, indicating that the system in the limit 
approaches the one with a NS tunnel junctionS^. On 
the other hand, in the S-lead dominant coupling regime 
shown in Fig|8l[b), the peak of zero bias conductance is 
split into two, which then move toward the opposite gap 
edges with keeping the unitary-limit value of 4e^//i when 
Ts increases. In this case, the SC-singlet is dominant at 
the QD in the equilibrium state, and the Andreev reso- 
nances, which originate from the exited doublet state for 
Fat = 0, emerge in the LDOS at the QD as shown in Fig. 
[5ja). The positions of the resonances are approximately 
given by iFg. We compare the voltage Vc that gives 
a peak in dl/dV with Fg in Fig. ^c). It is seen that 
eVc moves along the curve of F5 for large Fs/Fat, which 
confirms that the subgap peak in dl/dV results from the 
enhancement of the transport through the Andreev res- 
onances. 



A. Nonequilibrium transport for particle-hole 
symmetric case: ea/U — —0.5 

We now investigate the Coulomb interaction effects on 
the nonequilibrium differential conductance at a finite 
bias voltage. Since there are many relevant parameters 
in the system, we will divide our discussions into two 
cases. We first treat the simple case with a condition of 
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FIG. 9: (Color online) (a) Bias voltage dependence of the 
differential conductance for several values of f/: Ts/Tn ~ 
1, €d/U = -0.5, A/Fn = 0.5 and ksT/TN = 0.01. (b) 
Plots of the peak position, Va, in the dl/dV curve near the 
zero bias voltage and the renormalized coupling, Tn, which 
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FIG. 10: (Color online) (a) Local density of states at the QD 
for several values of V: Vs/Tn = 1, U/Fn = 20, €d/U = 
-0.5, A/Fjv = 0.5 and kgT/FN = 0.01. (b) Peak position of 
subgap resonances in (a) as a function of V. 



particle-hole symmetry, ed/U — —0.5. More generic cases 
with arbitrary conditions for ed and U will be discussed 
separately in the next subsection in comparison with the 
experiments. 



1. Coulomb interaction effects for Fs/Fjv — 1 

Let US start with a system with the symmetric cou- 
plings for tunneling, Ts/Tn = 1, which may help us to 
imagine what is essential in the nonequilibrium transport 
in the interacting QD. Figure IHJa) shows the differential 
conductance as a function of the bias voltage for several 
values of the Coulomb interaction, U. According to the 
analysis in the previous section (see Fig. [T]) , with increas- 
ing U, the system enters the N-lead dominant coupling 
regime where the Kondo singlet state becomes dominant. 
In this Kondo regime, several peaks appear at subgap 
voltages. We refer to the two sharp peaks near the zero 
bias voltages as Peak A and the two broad peaks at higher 
voltages as Peak B in Fig. ^a). Although the heights 
of the peaks are suppressed, both of Peak A and Peak B 
become prominent for large U. Note that the U depen- 
dence of the position of Peak A is different from that of 
Peak B; Peak A approaches the zero bias voltage with 
increasing U, whereas Peak B slightly shifts toward the 



gap edge. This fact implies that these two kinds of sub- 
gap peaks in dl/dV have different origins. We will show 
below that Peak A originates from the interplay between 
the Kondo effect and the Andreev reflection at a finite 
bias, while Peak B comes from the Andreev bound states 
at the QD. 

Let us focus on Peak A. Figure [HJb) shows the com- 
parison of the position of Peak A, denoted as Va, and the 
renormalized N-lead coupling Tn defined at V — T = 0. 
It is seen that the value of e Va approaches F n when the 
system enters the Kondo regime with increasing U. Since 
Fat is the characteristic energy scale of the Kondo ef- 
fect, which approximately gives the width of the Kondo 
resonance, Peak A is related to the Kondo effect. The 
emergence of the Kondo effect is seen in the bias voltage 
dependence of the LDOS at the QD shown in Fig. [TUl 
Although the LDOS for [/ = is not changed by the 
bias voltage, it is affected via the self-energy for finite U. 
In particular, the LDOS in the Kondo regime substan- 
tially changes its form under a finite bias voltage. Figure 
shows the LDOS at the QD for U/Tn = 20 and 
Ls/Lat — 1. For V — 0, there is a sharp Kondo resonance 
at the Fermi energy. With increasing V , the position of 
the Kondo resonance follows the chemical potential of the 
N-lead, fi^ — eV, suggesting that the Kondo screening 
of the local moment is mainly caused by the normal lead. 
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FIG. 11; (Color online) Temperature dependence of the differ- 
ential conductance for Fs/Fiv = 1, U/Tn ~ 20, ed/U — —0.5 
and A/Fjv = 0.5. 
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A noticeable change in the LDOS at finite bias voltage 
(Fig. [TUf a)) is the appearance of the additional resonance 
which is located at the counter position of the ordinary 
Kondo resonance; the ordinary Kondo resonance has a 
shoulder structure for eV j IS. = 0.2, which is changed into 
an additional resonance for eV j IS. = 0.4. This consider- 
ation naturally suggests that the additional resonance is 
caused by the Andreev reflection through the ordinary 
Kondo resonance (referred to as Kondo-assisted Andreev 
reflection); an electron which comes from the N-lead 
Fermi surface reaches the S-lead via the ordinary Kondo 
resonance, and then it is converted as a hole via the An- 
dreev reflection process. Since the electron has finite en- 
ergy measured from the S-lead Fermi surface, eV, the re- 
flected hole also has the same energy. This interpretation 
clarifies why the position of the additional resonance is 
located at the counter position of the Kondo resonance. 
Note that the additional resonance discussed here was 
previously realized by Sun et al^, but was not discussed 
in detail, in particular, about its physical relevance to 
the transport properties. We will address this issue with 
the use of the renormalized couplings, and demonstrate 
that it indeed provides a source of the marked change in 
nonequilibrium transport properties. 

The positions of the Kondo and additional resonances, 
a;+ and a;_ are shown in Fig. llOf b) as a function of the 
bias voltage. For small V, aj+ and a;_ follow the dotted 
lines which denote the position of the chemical poten- 
tial of N-lead and its counter position, ±fiN- Hence the 
crossover of the LDOS from the single peak to the dou- 
ble peaks occurs at eV ~ F^v where the distance between 
the two peaks is approximately given by their width. By 
comparing the results of the LDOS with dl/dV, we find 
that the crossover voltage in the LDOS approximately 
corresponds to the one giving Peak A in the differential 
conductance. Summarizing all these results, we conclude 
that Peak A in dl/dV originates from the Kondo-assisted 
Andreev reflection at a finite bias voltage. 



FIG. 12: (Color online) Differential conductance as a function 
of bias voltage for several values of Fg: U/Tm ~ 20, ed/U — 
-0.5, A/Fjv = 0.5 and kBT/TN = 0.01. 



In order to further confirm our interpretation for Peak 
A in dl/dV, we calculate the temperature dependence 
of dl/dV as shown in Fig. [TT] With increasing tem- 
perature, two peaks near zero bias voltage, which are 
classified as Peak A, are smeared and absorbed into the 
broad peaks of Peak B. The temperature dependence of 
Peak A supports that it is due to the Kondo-assisted An- 
dreev refiection. The characteristic temperature around 
which Peak A is smeared coincides approximately with 
Pat/Ftv — 0.086. We note that the equilibrium quan- 
tity of Fjv characterizes both the position and the T- 
dependence of Peak A in the nonlinear differential con- 
ductance. 

In contrast, we can see that Peak B in dl/dV is not 
directly related to the Kondo effect according to its T- 
dependence in Fig. [TTJ although the two peaks labeled 
as Peak B decrease their heights with increasing T, the 
broad peak structure still exists even at T/Tn = 0.25 
(higher than the Kondo temperature). We indeed find 
that Peak B is solely controlled by the Andreev reflec- 
tion, but not by the Kondo effect. Since the origin and 
the physical implications of Peak B are naturally seen 
in the asymmetric limit of Fs/Fjv 3> 1, we will discuss 
the physical properties systematically in asymmetric cou- 
plings below. 



2. Coulomb interaction effects for Fs/Fjv 7^ 1 

Here, we address how the asymmetry of couplings 
(Ts/Tn 7^ 1) alters the nonequilibrium transport prop- 
erties. We start with the differential conductance dl/dV 
shown in Fig. [T^] for several values of Ts in the case of 
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strong interaction U/Tm = 20. Here, we only show the 
resuhs either in positive or negative V since dl /AV is 
symmetric with respect to = for td/U = — 0.5. The 
crossover between the different regimes occurs around 
Ts/Fat ~ 2 where dl/dV has its maximum value 4e^//i 
at zero bias voltage. The profiles of dl/dV in the N-lead 
and S-lead dominant coupling regimes are shown in Fig. 
n^ a) and (b), respectively. 

In Fig. [T2lb). the d//dT^ curves show properties analo- 
gous to those in the noninteracting S-lead dominant cou- 
pling case shown in Fig. |Hl[b); the peak of dl/dV moves 
toward the gap edge with increasing F5. The peak val- 
ues, however, depend on the Coulomb interaction and 
become smaller than Ae'^/h for large U. The suppression 
of the peak is attributed to the inelastic scattering ow- 
ing to the Coulomb interaction. However, in the extreme 
limit of Fs ^ F^v, U, the superconducting correlation 
dominates the QD and thus the Coulomb interaction ef- 
fects are reduced, leading to the suppression of the in- 
elastic scattering. For Fg/Fjv = 7.5, therefore, the peak 
value increases again. For later discussions, we refer to 
these peaks as Peak C. On the other hand, the results 
of dl/dV in the N-lead dominant coupling regime(Fig. 
[T2r a)) are a little bit complicated; the peak at zero bias 
voltage splits with decreasing F5, forming a double-peak 
structure both in the positive and negative half of the 
subgap voltage regions. The origin of these two peaks is 
the same as discussed above, so that we denote them as 
Peak A and Peak B. With decreasing F5, the position 
of Peak A hardly changes, while that of Peak B shifts 
toward the gap edge. With further decreasing Fg, both 
of Peak A and Peak B reduce their heights. Therefore, 
for F5 — 7- 0, dl/dV in the gap is completely suppressed 
as is the case for [/ = 0. 

In order to elucidate the origin of the peak formation, 
we plot the F5 dependence of the positions of Peak A, B 
and C, which are labeled as Va, Vb and Vc, in Fig. [r^a). 
Here, we take U as the energy unit. For comparison, the 
renormalized couplings, Tn and F5, are also plotted. For 
^s/U < (>)0.1, Fat is larger (smaller) than Fg, namely, 
the system is in the N-lead (S-lead) dominant coupling 
regime. It is seen that Va and Vc approach the values of 
Tpf and F5 in the limit of Fg — > and 00, respectively. 
The crossover between these two limits appears around 
Ts/U — 0.1, where the peak is located at F = with 
the unitary limit value, Ae^ /h. At the crossover point, 
the width of the zero bias peak is simply scaled by F = 
Fat = F5 as shown in Fig. [T^T b). It is seen that the 
differential conductance calculated for several choices of 
A and F5 quickly decreases around \eV\ = F. 

The above properties in the conductance are clearly 
understood in terms of the Andreev bound states. The 
open diamonds in Fig. [T^ a) denote the energy Wf, of 
the Andreev bound states at the QD for Fat = T = 0, 
which is obtained with the NRG calculations^"—. As 
mentioned in the previous section, the system shows a 
transition between the magnetic doublet and SC singlet 
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FIG. 13: (Color online) (a) Semilog plot of the peak positions 
of dl /dV in the N-QD-S system, Va and Vb, as a function of 
Ts/U. The other parameters are Fiv/C/ = 0.05, ed/U = —0.5, 
A/[7 = 0.025 and ksT/TN = 0.01, which are the same as 
in Fig. 1121 For comparison, the renormalized couplings in 
the N-QD-S system with V = T = 0, Tn and fs, and the 
energy of the Andreev bound states in the QD-S system with 
T — 0, Ub, are also plotted, (b) Semilog plot of dl /dV at the 
crossover point against V/T for several choices of A and Fs- 
The other parameters are the same as in (a). The value of F 
is 0.0551, 0.124, 0.157, 0.187 and 0.213 for each curve from 
top to bottom. In the low voltage region with eV < F, dl /dV 
fall into a single curve. Inset shows the same data of dl /dV 
as a function of eV/U. 



states for Fat = T = 0. The transition point is eval- 
uated as r'g^ /U ~ 0.129 from the condition uJb = 0. 
For Ts < (>)F^^, the doublet (SC singlet) becomes the 
ground state and the Andreev bound states originating 
from the SC singlet (doublet) appear in the LDOS at the 
QD (see also the discussion in Fig. [7]). For finite Fat, 
the local moment of the doublet ground state is screened 
by the electrons in the N-lead. Therefore, the transition 
change into the crossover between the Kondo singlet and 
the SC singlet. It is clearly seen in Fig. [T^ a') that Vb 
is indeed related to the Andreev bound states since Vb 
approximately coincides with the energy of the Andreev 
bound states. This is also the case for Peak C, and the 
difference between them comes from whether the ground 
state is the Kondo singlet (Peak B) or SC singlet (Peak 
C). We therefore reveals the origin of Peak B and C; 
when the energy corresponding to the Andreev bound 
states is externally supplied by the applied bias voltage. 
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FIG. 14: (Color online) False color-scale representation of 
|(d4,dj-)| as a function of V and Fs for Tn/U — 0.05, ed/U = 
-0.5, A/Fat = 0.5 and ksT/rN = 0.01. The filled circles, 
triangles and squares indicate the Fg dependence of cVa/A, 
Vb/A and Vfc/A. 



the weight of the excited state at the QD is increased, re- 
sulting in the enhancement of the Andreev reflection. Al- 
though the nontrivial correspondence between the peak 
position of dl / dV and the energy of the Andreev bound 
states has already been discovered by Deacon et a P^i^° 
via the experimental studies, to our knowledge, this is the 
first numerical calculation which systematically clarifies 
the correspondence in both coupling regimes by taking 
into account the Kondo effect. 

The difference between Vb and Vc can be more clearly 
seen in the superconducting pairing correlation at the 
QD, Figure [U shows the false color-scale repre- 

sentation of as a function of Ts and V. The 

filled circles, triangles and squares on the representa- 
tion indicate the Ts dependence of eV^/A, eVe/A and 
eVc/A shown in Fig. [T31 In the N-lead dominant cou- 
pling regime, e.g., in the case of Ts/U — 0.05, the su- 
perconducting correlation at the QD is weak at F = 
since the Kondo singlet is dominant. With increasing V , 
I (did^) I shows the peak at the bias voltage where Peak B 
is located. This result clearly indicates the enhancement 
of the weight of the SC singlet state at the finite bias 
voltage. On the other hand, in the S-lead dominant cou- 
pling regime, e.g., in the Ts/U — 0.3 case, it is seen that 
|(cZ^d^)| is relatively large at V — 0, and monotonically 
decreases with increasing V. In particular, a rapid de- 
crease of occurs around V = Vc, implying that 
the weight of the magnetic doublet state increases in the 
system instead of the SC singlet state at V = Vc- All 
these features are consistent with the above interpreta- 
tion of the peaks in the conductance. 

It is instructive to consider the LDOS for finite V to see 
the nature of the bound states. Since the V dependence 
of the LDOS in the N-lead dominant coupling regime has 
already been discussed in Fig. [TUl here we focus on the 
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FIG. 15: (Color online) (a) Local density of states at the 
QD for several values of V for Ts/U = 0.15, Tn/U = 0.05, 
ea/U = -0.5, A/Fjv = 0.025 and ksT/TN = 0.01. 



S-lead dominant coupling regime. Figure fTSl a) shows the 
LDOS at the QD for Ts/U = 0.15. In this S-lead dom- 
inant coupling case, there are the Andreev resonances 
corresponding to the excited doublet for Tn = 0. The 
double peaks of the LDOS in the equilibrium state are 
located at w/A ~ ±fs ~ ±0.5. For eV/A = 0.2, the 
resonances show little change as if the Coulomb interac- 
tion is absent in the system. For eV/A > 0.4, however, 
the bias voltage increases the distance between the two 
resonances and smears them. In order to see the V de- 
pendence of the resonances in detail, we plot the peak po- 
sition of the resonances as a function of V in Fig. [TSlb). 
Here, lu+ and oj- denote the positions of the peaks for 
positive and negative w, respectively. With increasing 
V, they move toward the opposite gap edges. In particu- 
lar, the peaks become sensitive to the change of the bias 
voltage around eV/ A ~ 0.3, which approximately cor- 
responds to the value of eVc/A. This feature indicates 
that the peaks tend to follow the chemical potential of 
N-lead and its counter position, ^jv and —^J.N, which can 
be regarded as a kind of pinning effect of the Andreev 
resonances. Note that the profiles of cj^ and — a;_ ap- 
proximately coincide with each other, implying that the 
resonances keep their symmetric structure with respect 
to the Fermi level of S-lead. Besides, in the N-lead dom- 
inant coupling regime, the pinning of the resonance be- 
comes more prominent as discussed in Fig. 1101 The origin 
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FIG. 16: (Color online) False color-scale representation of 
duj+/dfiN as functions of V and Fs for Fn/U = 0.05, ea/U — 
-0.5, A/U = 0.025 and kBT/TN = 0.01. The filled circles, 
triangles and squares indicate cVa/A, Vb/A and Vc/A. 



of the pinning is attributed to the Kondo effect, and is 
understood as foUows: in the S-lead dominant coupling 
case, the electron correlation is practically negligible in 
the low energy and low voltage region, except for the 
renormalization effects, since the SC singlet is dominant 
at the QD. With increasing V, however, the weight of 
the magnetic doublet state increases in the system near 
V ^ Vc- Then the resulting doublet state is screened 
by the electrons in the N-lead owing to the Kondo effect, 
leading to the pinning of the resonances. 

In order to further investigate the pinning of the res- 
onances in the LDOS, we show the false color-scale rep- 
resentation of dw+/d/LtAr in Fig [161 Since dw+/d/XAr be- 
comes large if the resonance of the LDOS in the pos- 
itive Ld region follows the chemical potential of the N- 
lead, its value gives an estimate of how strong the Kondo 
correlation is. In the N-lcad dominant coupling regime, 
Ts/C^ ^ 0.1, the Kondo pinning effect is suppressed with 
increasing V since the bias voltage destroys the Kondo 
singlet state. In particular, da;+ /d^N is rapidly decreases 
when the value of the bias voltage approaches Vb be- 
cause the weight of the SC singlet is increased at the 
voltage. In contrast, in the S-lead dominant coupling 
regime, du!+/d^N increases with increasing V from zero 
and takes a peak at a finite bias voltage where Peak 
C appears in dl/dV. Therefore, it is intuitively under- 
stood that the Kondo correlation is enhanced by the bias 
voltage around V = Vc- With further increase in V, 
the Kondo correlation is weakened again by the applied 
bias voltage, which in turn leads to the suppression of 
du!+/d^N- The different features in da;+/d^Ar atV = Vb 
and V — Vc reflect the difference in the origin of the An- 
dreev bound states. Therefore, all the results of dl/dV, 
\{d^d/\-)\ and duj^/dfiN are consistent with the scenario 
that the weight of the excited states is enhanced when 




FIG. 17; (Color online) False color-scale representation of 
the magnetization at the QD as functions of V and Fs for 
Fn/U = 0.05, ea/U = -0.5, A/U = 0.025 and ksT/FN = 
0.01. 

the strength of the bias voltage coincides with the energy 
of the excited states. 



3. Comparison with mean-field results 

In order to clarify the role of electron correlations at a 
finite bias voltage, it is instructive to compare the present 
results with the mean-field approximation in U. In the 
mean-field approximation, there are two types of the so- 
lutions: the magnetic and nonmagnetic ones. The mag- 
netic solution appears only in the strong Coulomb in- 
teraction case, as known in the Anderson model2^. Al- 
though the magnetic solution does not describe the cor- 
rect physics at zero temperature, we can infer the en- 
hancement of the magnetic correlations at the QD via its 
existence. Hence, the mean-field analysis highlights the 
importance of the correlation effects, as discussed in the 
Anderson model out of equilibrium^^ . 

Figure [T71 shows the local magnetization m = ((%) — 
(n^))/2 at the QD as a function of V and Fg. In the 
case of 0.25 < Ts/U < 0.3, it is noteworthy that the 
magnetic solution exists only at a finite bias voltage; the 
bias voltage induces the local moment at the QD. With 
further increase in V , however, the magnetic solution dis- 
appears again. Such a reentrant behavior is not found in 
the nonequilibrium N-QD-N system for ed/U = — 0.5^i. 
Accordingly, we attribute the reentrant behavior to the 
competition/cooperation of the magnetic and supercon- 
ducting correlations at a finite bias voltage. The reen- 
trant behavior of the magnetic boundary clearly indicates 
that the magnetic correlation is enhanced by the bias 
voltage. Of course, the magnetic state is an artifact of 
the approximation and should be replaced by the Kondo 
singlet state, but the reentrant behavior of the boundary 
is consistent with the enhancement of the Kondo corre- 
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Andreev bound states, except around the boundary line. 
At the boundary, the conductance changes in a discon- 
tinuous manner, implying that the physics related to the 
Andreev bound states can be qualitatively understood 
without quantum fluctuations^^. 

On the other hand, quantum fluctuations induce the 
intriguing phenomena a la Kondo, which cannot be de- 
scribed at the mean-field level in the region where the 
magnetic and superconducting correlations compete with 
each other. Figure [rsFb) shows the MPT results of 
dl/dV. In addition to the good correspondence between 
the ridge of dl/dV and the energy of the Andreev bound 
states, the anomalous enhancement occurs around the 
transition point of the ground state in the QD-S system, 
which is caused by the Kondo-assisted Andreev reflec- 
tion. The quantum fluctuation effects are significant at a 
finite bias voltage in the vicinity of the transition point; 
The Kondo-assisted Andreev reflection, as well as the 
pinning effect of the Andreev resonances discussed in the 
previous section, is seen only at a finite bias voltage and 
smeared away from the transition point. Applying the 
bias voltage resembles increasing temperature concern- 
ing the enhanced weight of the excited states. 

B. Nonequeliburium transport for generic (e^, V) 
cases: comparison with experiments 

In the remainder of the section, we discuss the differ- 
ential conductance as functions of the energy level at the 
QD, ed, and the bias voltage, V. Note that can be eas- 
ily controlled by the gate voltage in actual experiments. 
We compare the analysis done here with the recent ex- 
periments qualitatively in good agreement 



FIG. 18: (Color online) (a) Mean-field approximation results 
of the differential conductance: Tm/U = 0.05, ed/U = —0.5, 
A/U = 0.025 and fesT/Fiv = 0.01. The energy of the An- 
dreev bound states, W(,/A, as a function of ea/U is denoted 
by diamonds. The black line indicates the boundary of the 
magnetic solution, (b) MPT results of the differential con- 
ductance. The parameters are same as those in (a). 

lation deduced from the MPT calculation. 

We also calculate the differential conductance, dl/dV, 
with the mean-field approximation. Figure llSf a) shows 
the false-color scale representation of dl/dV as functions 
of eV/A and the logarithm of Ts/U. The black solid 
line indicates the magnetic boundary, and the diamond 
denotes the energy of the Andreev bound states. In the 
region where the magnetic solution exists, the ridge posi- 
tion of dl/dV features a curve similar to the energy of the 
Andreev bound states. On the other hand, in the large 
Ts region where the non-magnetic solution only exists, 
the ridge is located around the gap edge and coincides 
with the Andreev bound states. As a result, the mean- 
field solutions still capture the essential nature of the 



1. Transport via Andreev bound states 

Let us first focus on the simple cases where the Kondo- 
assisted Andreev transport is not so important, for which 
we can highlight the crossover behavior in the transport 
due to the Andreev resonances. 

Figure [19] is the false color-scale plot of dl/dV as a 
function of ed for several values of Fg. We also show the 
ed dependence of the energy of Andreev bound states, 
which is calculated for Fat = T = with using the NRG 
methodSS. In Fig. fTWa) and (b), the system is in the S- 
lead dominant coupling regime. It is clearly seen that the 
value of the bias voltage where dl/dV takes a peak value 
coincides approximately with the energy of the Andreev 
bound states. The peak ridges correspond to Peak C 
defined in the previous subsection. For Ts/U = 0.16, the 
system shows a crossover around ed/U = —0.5 and the 
conductance has the maximum value, Ae^/h (center of 
the figurefTWc)). Note that the energy of Andreev bound 
states does not touch the zero axis since Ts/U = 0.16 
characterizing the crossover is slightly different from the 
exact transition point. However, the peak positions still 
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FIG. 19; (Color online) False color-scale representation of the differential conductance as functions of V and for Tn/U — 0.1, 
A/U = 0.05 and ksT/VN = 0.01: (a) Ts/U = 0.5, (b) Vs/U = 0.25, (c) Ts/U = 0.16, (d) Ts/U = 0.1, (e) Ts/U = 0.05. The 
energies of the Andreev bound states, tot/ A, are also plotted. 



approximately correspond to the energy of the Andreev 
bound states. 

With decreasing Fg, the system enters the N-lead dom- 
inant coupling regime as shown in Fig. 119( d). In the 
case of Fat — T = with the other parameters same as 
in (d), the magnetic doublet becomes the ground state 
around ed/U = —0.5, and a transition occurs, where 
the softening of the Andreev bound states occurs around 
ed/U ~ —0.5 ± 0.32. After the transition, the Andreev 
bound states move toward the gap edges. For finite F^r, 
the magnetic doublet is changed to the Kondo singlet, 
which is realized around the center of the figure. Re- 
garding dl/dV, the peaks corresponding to the Andreev 
bound states are not so clearly seen as in the S-lead dom- 
inant coupling regime. This is because Fjv is not small in 
comparison with F5. The peaks around the symmetric 
point are especially indistinct because there are two kind 
of peaks of dl/dV: Peak A originating from the Kondo- 
assisted Andreev reflection and Peak B corresponding to 
the Andreev bound states. With further decreasing F5, 
the bare value of the coupling of the S-lcad, Fg, becomes 
smaller than that of the N-lead, Fat. For F^r > F5, these 
peaks almost disappear in the subgap voltage. However, 
there is a remnant of the Kondo-assisted Andreev reflec- 
tion, which makes a dip at y = in Fig. [TWe). For 
Fg — > 0, the remnant is also diminished and there are 
only peaks at the gap edges as in the noninteracting N- 
lead dominant coupling regime. 

As a result, both for Tn/Ts and 00, the subgap 



conductance, except at the gap edges, is completely sup- 
pressed for any values of e^. A wide variety of patterns of 
the differential conductance result from the competition 
between the Coulomb interaction and the superconduct- 
ing correlations at the QD. 

It is to be noted here that the pronounced gap-edge 
peaks in the cases of Fat Fg and F^r ^ Fg, and the 
prominent peaks corresponding to the Andreev bound 
states in the S-lead dominant regime are qualitatively in 
agreement with the recent experimental resultsSSJS,. In 
the experimen t^^'^" , the softening of the Andreev bound 
states in the N-lead dominant coupling regime is also 
observed in the distinctive peaks in the dl /dV measure- 
ment. However, we do not find such a distinctly sep- 
arated peak in the case Fjv > Pg. We will discuss the 
visibility of the peaks separately below. 

We also calculate the (e^, V) dependence of as 
shown in Fig. [20l For Ts/U = 0.5, it is seen that \ {did^)\ 
is large around the center of the figure and decreases 
with increasing since the dominant SC-singlet state, 
which consists of the superposition of doubly-occupied 
and empty states, is weakened in the empty or doubly- 
occupied region. For finite V, |(c?^d^)| decreases rapidly 
around the voltage corresponding to the energy of the 
Andreev bound states of a doublet character. For small 
Fs, the peak of |(c?^d^)| is divided into two which are 
located around the voltage corresponding to the energy 
of the Andreev bound states of a SC-singlet character. 
In spite that the electron and hole components of An- 
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FIG. 20: (Color online) False color-scale representation of |{d4_d-|-)| as functions of V and for Fjsi/U = 0.1, A/U = 0.05 and 
kBT/VN = 0.01: (a) Vs/U = 0.5, (b) Vs/U = 0.25, (c) Vs/U = 0.16, (d) Vs/U = 0.1, (e) Vs/U = 0.05. The energies of the 
Andreev bound states, aj(,/A, as a function of ed/U are also plotted. 



dreev bound states have the same energy, there is asym- 
metry in dl/dV and as a function of V. For 
instance, \{did^)\ is strongly suppressed in the right top, 
in comparison with the one in the right bottom in Fig. 
[501 Moreover, the asymmetric feature becomes promi- 
nent with decreasing Fg. In order to explain how the 
asymmetry emerges, let us consider the F5 — limit. 
Figure [21] shows the schematic phase diagram of the N- 
QD system with strong U. The Kondo correlation is 
dominant for ed/U — —0.5 and ^ = 0, and is weakened 
with increasing je^ -I- 0.5|. Note that V just shifts the 
chemical potential of the N-lead, so that only one of the 
two parameters, eV and e^, becomes relevant; the sys- 
tem stays in the same state along fij^ = line in the 
figure. For —0.5 -I- A/U < ed/U, therefore, the Kondo 
correlation is enhanced if we fix and increase V. This 
enhancement of the Kondo correlation would occur in the 
N-QD-S system with small F5, giving rise to the asym- 
metric patterns. 



2. How to observe Kondo- assisted Andreev transport 

The subgap peak structure in the differential conduc- 
tance in the N-lead dominant coupling regime is a bit 
more complicated than that in the S-lead dominant cou- 
pling regime because both of the Kondo-assisted Andreev 
reflection and the Andreev bound states could affect the 
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FIG. 21: (Color online) Schematic phase diagram of the N- 
QD-S system for Fs = and large U. 



conductance profile as seen in Fig. 1191 Therefore how we 
can observe the conductance peaks in the N-lead domi- 
nant coupling regime depends sensitively on the ratio of 
the bare coupling strengths. We find that if the N-lead 
dominant system is in the condition Tpj/Ts < 1, that is 
if Fat/Fs > 1 and F^r/Fg < 1, the conductance peaks 
become prominent in the subgap voltage. We show some 
results below in the case satisfying this specific condition. 

Figure[221sliows the differential conductance for several 
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FIG. 22: (Color online) Differential conductance as a function 
of V for several values of (Fs, A): Tn/U = 0.05, ed/U = 
-0.5, and ksT/VN = 0.01. 



values of (F, A) and ed/U — —0.5. Note that the system 
is in the N-lead dominant coupling regime for any choices 
of (F, A) (see also the phase diagram of Fig. [7]). The 
two peaks in the vicinity of ^ = are related to the 
Kondo-assisted Andreev reflection and the other two at 
eV/ Is. ~ ±0.5 originate from the Andreev bound states. 
These four peaks become sharper and more prominent 
from the top to bottom lines. In particular, the Kondo- 
like peaks approach the zero bias voltage since the peaks 
are located at eV ~ ±Tm and Fjv decreases from top to 
bottom. Therefore, we conclude that for [/ ^ F5, A ^ 
F^v, only a single peak around F = 0, instead of the 
two peaks at finite bias voltages, may be observed in real 
experiments. 

We next look at the dependence of the N-QD-S sys- 
tem for Tm/U ^ 0.05, TsjU = and A/C/ = 0.075 
(Fig. [23]). In this case, Fjv > is satisfied at the 
particle-hole symmetric point {ed/U — —0.5) in addi- 
tion to the condition of Fjv/Fs < 1, which leads to the 
crossover in the dominant couplings may occur away from 
the symmetric point as shown Fig. ISBT a) . The crossover 
points are approximately coincide with the doublet- 
singlet transition points denoted by arrows for F^r = 0. 
In the particle- hole symmetric case, Tm/^n — 0.1 char- 
acterizes the position and the temperature dependence of 
the Kondo-type peaks in dl/dV. For ksT < Fjv, both of 
the Kondo-assisted Andreev reflection and the Andreev 
bound states at the QD contribute to the nonlinear elec- 
tron transport. Figure [23l b) shows the Sd dependence of 
dl/dV at /cbT/Fat = 0.01. Regarding the central two 
peaks related to the Kondo effect, one of the peaks in- 
creases its height away from ed/U = —0.5, whereas the 
other disappears. We note that the curve of dl/dV ends 
at certain voltages for €d/U — —0.4, 0.3, 0.2 and —0.1 
since we cannot get the convergent solution in the frame- 
work of MPT. Therefore, the movement of the one of 
the two peaks owing to the Andreev bound states is not 
clear in this figure, unfortunately. Nevertheless, we can 



see that the other peak moves to lower bias voltages and 
merges into the single peak with the prominent peak re- 
lated to the Kondo effect. With further increasing ed, the 
merged peak moves toward the gap edge and there is no 
peak in the subgap voltage. 

A quantitative comparison of the peak positions with 
the energy of the Andreev bound state, Wb/A, is shown 
in Fig. [23jc). Softening of the Andreev bound states 
clearly emerges in the peak structure of the differen- 
tial conductance. The softening reflects the transition 
of the ground states in the case that Fat = 0. For 
—0.75 <; €d/U < —0.25, there are two Kondo-assisted 
Andreev reflection ridges around V = 0. The two 
ridges are separated around ed/U — 0.5. In the case of 
^d/U ^ —0.75 or —0.25 > Cd/U, the Kondo-type ridges 
are not observed since the system is away from the Kondo 
regime, and there appear only the gap-edge peaks related 
to the Andreev bound states. The dashed lines in the fig- 
ure denote the conditions where the level of the QD and 
the chemical potential of the N-lead satisfy the equations; 
^J'N = + U and fJ,N = ^d- In spite that the Andreev 
bound states have the symmetric energy spectrum, the 
large asymmetry is found in dl/dV as a function of V, 
except for €d/U 7^ —0.5. The asymmetry reflects the fact 
that the Kondo correlations appear differently depending 
on the sign of the bias voltage, as discussed in Fig. [21] 

It is remarkable that the overall features of the differ- 
ential conductance are consistent with those in the re- 
cent experiment"^*^, including the observation of the fln- 
gerprints of two kinds of peaks; Kondo-type peaks at 
eV/A ~ and the peaks due to the Andreev bound 
states. There still seems to be a small discrepancy be- 
tween the theory and the experiment^''. In our theory, 
the Kondo-type ridges are separated at the particle-hole 
symmetric point, while a single Kondo ridge is observed, 
instead of the two ridges, in the experiment. We believe 
that the single ridge is a consequence of the special condi- 
tion, U 3> r5,A ^ F^v, used in the experiment. In this 
case, the two Kondo-type peaks for ed/U = —0.5 would 
be located near = and overlap with each other, as 
discussed in Fig. [22] Hence, the two Kondo ridges could 
be also connected at the particle-hole symmetric point in 
such a special condition. 

Finally, some comment are in order for the tempera- 
ture dependence. In Fig. [23l d). we present dl/dV at a 
higher temperature ksTx: /T n = 0.1, which is compara- 
ble to F^v. The Kondo-type ridges are smeared with in- 
creasing temperature, but the peaks due to the Andreev 
bound states show little change. Therefore, at higher 
temperatures, there would be only the crossover behav- 
ior of the fingerprint of the Andreev bound states in the 
dl/dV measurements. 



V. SUMMARY 

In this paper, we have theoretically investigated the 
nonequilibrium electron transport through a quantum 
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FIG. 23: (Color online) (a) Renormalized coupling strengths as a function of td for Fn/U = 0.05, Ts/U = 0.1, A/U = 0.075 
and T = V = 0. (b) Differential conductance as a function of the bias voltage for several values of ed/U and ksT/TN = 0.01. 
From bottom to top, ed/U = —0.5, —0.4, —0.3, —0.2, —0.1, 0, 0.1, 0.2 and 0.3. The other parameters are the same as in (a), 
(c) False color-scale representation of the differential conductance as a function of V and ed- The parameters are the same as 
in (b). The energies of the Andreev bound states, uJb/A, are denoted as diamonds. The bound states exist for Fjv — T = 0. 
The dashed lines indicate the resonant conditions; pjv = and /ijv = ed + U. (d) At ksT /Tn = 0.1 which is nearly equal to 
that of Fjv /Fjv, the central Kondo-enhanced Andreev ridges are suppressed. 



dot coupled to the normal and superconducting leads 
with particular emphasis on the interplay between the 
Kondo effect and the superconducting correlations. For 
this purpose, we have developed the modified second or- 
der perturbation theory in Keldysh-Nambu formalism 
under nonequilibrium steady-state conditions. We have 
confirmed that this method is indeed efficient for analyz- 
ing the nonequilibrium electron transport in the present 
system. 

It has been shown that the renormalized couplings be- 
tween the leads and the dot in the equilibrium states are 
the key quantities that correctly describe nonequilibrium 
transport properties. In particular, the enhancement of 
the Andreev transport occurs via a Kondo resonance at 
a finite bias voltage, giving rise to an anomalous peak 
structure in the differential conductance, whose position 
is determined solely by the above-mentioned renormal- 
ized parameters. This peak formation is a remarkable 
example of phenomena that are indeed caused by the 
interplay between the Kondo and superconducting cor- 
relations. A pinning effect of the Andreev resonances to 
the Fermi level of the normal lead and its counter level 
also evidences the interplay of the above two types of 
correlations. Moreover, it has been shown that the en- 
ergy levels of the Andreev bound states give rise to an 
additional peak structure in the differential conductance 
in the strongly correlated N-QD-S system. 

We have demonstrated that the above characteristic 



features of nonequilibrium differential conductance ob- 
tained from our calculation are qualitatively in agreement 
with those observed in the recent experiment o^^i'^° . Fi- 
nally we note that there still exists a small discrepancy 
between the theory and the experiments; the Kondo-type 
ridges are separated in our theory (Fig. I23p . while a sin- 
gle Kondo ridge is observed experimentally.^^ We believe 
that the single ridge is a consequence of a specific con- 
dition employed in the experiments, U S> F5,A 3> Fjv 
and that the clearly-separated Kondo ridges predicted in 
this paper will be observed if the proper conditions for 
the system parameters are prepared experimentally. 
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